Fermion Tunneling from Dynamical Horizons 
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Abstract 



The instability against emission of fermionic particles by the trapping horizon of an 
evolving black hole is analyzed using the Hamilton- Jacobi tunneling method. This method 
automatically selects one special expression for the surface gravity of a changing horizon. 
The results also apply to point masses embedded in an expanding universe. As a bonus of 
the tunneling method, we gain the insight that the surface gravity still defines a temperature 
parameter as long as the evolution is sufficiently slow that the black hole pass through a 
sequence of quasi-equilibrium states, and that black holes should be semi-classically unstable 
^ c"| even in a hypothetical world without bosonic fields. 

PACS numbers: 04. 70.-3,04. 70.Bw,04. 70.Dy 
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m . 1 Introduction 

In a previous paper [1] we considered the quantum instability of dynamical black holes using 
a variant of the tunneling method introduced by Parikh and Wilczek [2] , according to which 
the Hawking effect can be considered as a kind of tunneling transition through the horizon 
of the black hole (BH). Probably, for such matters it was first applied in condensed matter 
physics, specifically to exotic phenomena involving the analogue of horizons in superfluid 
3 He films [3]. The method was refined and extended to more general cases in [4] and others 
papers as well [5] , including the back reaction effects [6] and its extension to non-commutative 
Schwarzschild space-time [7] . For criticism and counter criticism see also [8] , and [9] for the 
relation between tunneling and thermodynamics. 

The Hawking's effect should not be confused with another tunneling process, the escape of 
the particle to infinity through the potential barrier surrounding the horizon, and about 
which the Parikh- Wilczek method has nothing to say. 

The variant mentioned above is the Hamilton- Jacobi method introduced in [10,11], so called 
after the appearance of a complete comparison analysis with the Parikh- Wilczek method done 
in [12]. The tunneling method provides not only new physical insight to an understanding 
of the black hole radiation, but is also a powerful way to compute the surface gravity for a 
vast range of solutions. Now for non stationary black holes, termed dynamical black holes 
in [1], things are not so simple and even the possibility of Hawking radiation is in principle 
questionable, since in general a changing horizon is not a null hypersurface, although it is still 
one of infinite red shift. In particular, several definitions of the surface gravity for evolving 
horizons have been proposed in the past, all fitting the first law of BH mechanics more or less 
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equally well. A comparison is discussed throughly in [13]. Among these the one proposed 
implicitly by Visser [10] and more explicitly by Hayward and Kodama [14, 15] is the most 
interesting to us, as it is the one the tunneling method leads to. Still different, otherwise 
reasonable, resuls are advocated for expanding cosmological black holes [16,17], but we shall 
not dwell about them here. 

As we said above, a way to understand Hawking radiation is by means of tunneling of 
particles through black-hole horizons. Such tunneling approach uses the fact that the WKB 
approximation of the tunneling probability for the classical forbidden trajectory from inside 
to outside the horizon is: 

roce-« Im/ , (1.1) 

where / is the classical action of the trajectory, to leading order in h. What kind of particles 
do we expect to find in the Hawking radiation spectrum of a black-hole? In principle, 
all the Standard Model particles. However, most of the calculations in literature have been 
performed just for scalar fields, except in [18], where a detailed study of spin one-half emission 
was considered for stationary black holes, and [19] for the special case of the BTZ black hole. 
What we are going to check is that the tunneling approach via the Hamilton- Jacobi method 
is consistent with the Kodama-Hayward prediction even for 1/2— spin particles in dynamical 
black holes and that, as a consequence, they should be semi-classically unstable even in a 
world with no bosonic fields. 



2 Bardeen-Vaidya space-times and cosmological black 
holes 

Let us consider the Bardeen - Vaidya metric (BV) in D = 4 dimensions [20, 21] 

ds 2 = ~e 2 ^ r - v) A(r, v)dv 2 + 2e^ r ' v Uvdr + r 2 dcj 2 , (2.1) 

where v is an advanced time null coordinate, dto 2 = d6 2 + sin 2 8 d<p 2 and, in the simplest case 
we have in mind, A(r,v) is just 

A(r,v) = l-^. (2.2) 
The inverse metric is given by: 
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(2.3) 



This class of metrics possesses - under very general conditions - a trapping horizon (TH), 
as defined by Hayward [14], which in the present case is given by the equation A(r, v) = 0. 
This defines a curve r — ru(v) giving the location of the apparent horizon; the quantity 
E = m(rH{v),v) is the Misner- Sharp mass [22] of the horizon, in term of which the horizon 
will be trapping if m {th,v) < 1/2, a prime denoting the radial derivative. The TH proved 
to be the key concept in such matters as particles emission and gravitational entropy (for a 
very recent discussion see [23]), and can be applied successfully even to higher dimensional 
Vaidya space-times [24]. 

According to the Kodama-Hayward theory, to such TH is associated a geometrical surface 
gravity 

^) = J__m'(r H ,v) 

We see the meaning of the trapping condition: it ensures the positivity of the surface gravity. 
The second example we are interested in is the McVittie solution [25] for a point mass in 
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a Friedmann-Robertson- Walker flat cosmology. In isotropic spatial coordinates it is given 
by [25] 

ds 2 = -A(p, t)dt 2 + B(p, t) (dp 2 + p 2 duj 2 ) (2.5) 

with 



A(p,t) = 



1 



a{t)p 



1 



1 - 



a(t)p 



B(p,t) = a(ty 

a(t)p _ 

When the mass parameter m = 0, it reduces to a spatially flat FRW solution with scale factor 
a(t); when a(t) — 1 it reduces to the Schwarzschild metric with mass m. In four dimensions 
this solution has had a strong impact on the general problem of matching the Schwarzschild 
solution with cosmology, a problem faced also by Einstein and Dirac. Besides McVittie, it 
has been extensively studied by Nolan in a series of papers [26]. To put the metric in the 
general form of Kodama theory, we use what may be called the Nolan gauge, in which the 
metric reads 

ds 2 = - (A 8 - H 2 {t)r 2 ) dt 2 + A^dr 2 - 2A~ 1/2 H(t)r drdt + r 2 dw 2 (2.6) 

where H(t) — a/a is the Hubble parameter and, for example, in the charged case, A s — 
1 — 2m /r + q 2 jr 2 . In passing to the Nolan gauge a choice of sign in the cross term drdt 
has been done, corresponding to an expanding universe; the transformation H(t) — > —H(t) 
changes this into a contracting one. In the following we shall consider q = 0; then the 
Einstein-Friedmann equations read 

3H 2 = 8tt P , 2A- 1/2 H{t) + 3H 2 = -8np . (2.7) 

It follows that A s = 0, or r = 2m, is a curvature singularity similar to r = in FRW models, 
namely it is a big bang singularity. When H — one has the Schwarzschild solution. The 
term H 2 r 2 in the metric strongly resembles a varying cosmological constant; in fact if H is 
constant the metric reduces to the Schwarzschild-de Sitter solution in Painleve coordinates. 
The trapping horizon is the root of A s = H 2 rj^, and is time dependent. The Misner-Sharp 
mass and the geometrical surface gravity are, respectively, 

E = m + ±H(t) 2 r 3 H , (2.8) 

(2.9) 

Note that E — r#/2. In the static cases everything agrees with the standard results. The 
surface gravity has an interesting expression in terms of the sources of Einstein equations 
and the Misner-Sharp mass. Let T^d be the reduced trace of the stress tensor in the space 
normal to the sphere of symmetry, evaluated on the horizon. For the Vaidya-Bardeen metric 
it is, by Einstein's equations, 

i civn 

T 2D = TZ+T r r = -- — (2.10) 

Inrn Or \r=r H 

For the McVittie's solution, this time by Fredmann's equations (|2.7p one has 

T„ = - P+P = -±(3H> + JL^ (2.11) 
We have then 

k = + 2nr H T 2D . (2.12) 

It is worth mentioning the pure FRW case, i.e. A s = 1, for which n(t) = —H(t) — H/2H. We 
feel that these expressions for the surface gravity are non trivial and display deep connections 
with the emission process. We now proceed to analyze the Dirac equation in BV space-time. 
Since the Hawking's effect is a phase phenomenon, which is why the tunneling method works 
so well, we expect the spinor amplitudes to play a minor role. We verify this for the BV 
"shining star" , but use a more straightforward method for the McVittie's solution. 
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3 Dirac equation for BV metric 

The Dirac equation is, 

(VI> M + ^)v(«,r,e,0)=O, (3.1) 

where 

D„ = d » + \K ^ ( 3 - 2 ) 

S Q/ 3 = -[7a, 73]-. (3-3) 

The 7*'— matrices satisfy the Clifford algebra, 

[ la , lfj ] + = 2g aP l, (3.4) 

where I is the (4 x 4)-identity matrix. In order to get the Dirac-7 matrices for the BV metric 
at hand, we define firstly a tetrad of orthogonal vectors s.t. 

VabE^Et = 9tlv . (3.5) 

The convention is that (first) latin indices are Minkovskian so they run over (0, 1, 2, 3); greek 
indices are coordinate indices so they run over (v, r, 9, 0). 

Of course there are many different tetrads, but the simplest choice is the following: 

E a v = (eVA, 0,0,0), (3.6) 

K = (-vTvT°'°)' (37) 

E% = (0,0,r,0), (3.8) 
El = (O,O,O,rsin0). (3.9) 
The 7 M — matrices are expressed in terms of the tetrad in the following way: 

l^laE^. (3.10) 

With our choice (|3.6|) - (|3.9p . it turns out that 

(3.11) 
(3.12) 

(3.13) 
(3.14) 



lv 




lr 


= 7i (7l ~ 


16 




l<t> 


= r sin 9 73 ; 



l v = -=(7i-7o), (3.15) 



and 



7 r - VA 7l , (3.16) 

1° = -72, (3.17) 

r 

I 4, = ^^73- (3.18) 
We can also express the 

7 5 d M Y 1° 1* = -^-— D (I - 7o7i) 7273- (3.19) 
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All what we need now is a convenient representation of the Dirac-7 a matrices satisfying 
[7a,7h]+ = 2r] ab with 77 = diag(-l, 1, 1, 1). 



7o 
7i 

72 

73 
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-i 

a 3 

cr 3 

o-i 

01 

<7 2 

(72 



The a- matrices are the Pauli matrices satisfying the usual relations: 



aw = l( 2x %- + 
By virtue of (|2.ip and 



result in the footnote below^ for the sake of completness. Next, we should calculate the S Q /j 
defined in (|3.3|) : a task that Maple can do very quickly; in the end, it remains to evaluate 
the combination of such results, namely 



i,j,k = 1,2,3. 
one calculates the T a ^ symbols required in 



(3.20) 
(3.21) 
(3.22) 
(3.23) 

(3.24) 
We list the 



-V = 0. 



(3.25) 



Let us employ the following ansatz for the spin-up Dirac fielqj: 

¥(u,r,#,0) T = 





to(v,r,6,<t>) 




exp 



-J T (u,r,0,. 



(3.26) 



Plugging the ansatz (|3.26p into Dirac equation (|3.25[) . it turns out that the term in square 
brackets is of order 0(K). Thus, we do not need to work out its precise form, since in the h — > 
limit it vanishes. To leading order in h equation (|3.25|) becomes = (h$ + to) + 0(h), or 
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-_If(v,r,6,i 








(3.27) 
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2 We shall perform a detailed analysis only for the spin-up case, being confident that mutatis mutandis every- 
thing applies in the same way to the spin-down case. 
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Thus, we get the following equations: 



(e _ *0„) J T - -= (e-^fllu) /t - V^^t + iEm = 0; 
^4 y A 



-[del, 
r \ sin ( 



d v )h 



<•/,/• l=D 

in 



■- a e / T + — -o^if 

r \ smO 



e'^dv) 7 T - VAEd r Ii + iflm = 0; 



0, 



(3.28) 
(3.29) 
(3.30) 
(3.31) 



and the Kodama vector K = e~^d v has been put in evidence throughout. K plays the role 
of the Killing vector <9 t for dynamical black-holes. Therefore, it makes sense the following 
ansatz for the action: 



'1 



dve^ r ' v) E + W(r) + J{6,4>), 



which inserted into (|3~28f - ([3~3"Tj) gives 
v : 



-=(iE - Q)E - \/MW'(r) + iEm = 0; 
"A 



r \ sine> 



-=(E + ifL)E - \fAEW'(r) + ittm = 0; 
^4 



M0,<l>) 



-,J<b 



0. 



(3.32) 

(3.33) 
(3.34) 
(3.35) 
(3.36) 



(|3.34p and (|3.36|) imply that J(9, cf>) is a complex function. The same solution for J is 
obtained for the spin-down case, then its contribution to the rate emission T cancels out and 
we can forget about it. As regard the remaining equations, we have 

1. E = zO, then 

2E 



W'(r) 



(3.37) 



A(r,v)' 
2. E = — ifl, then 

W'(r) = 0, (3.38) 

both in the massless and massive cases. Solution ()3.38|) has not to be considered surprising, 
since 



7 5 ¥ = 





5 + iil 




(3.39) 



meaning that case 2. corresponds to the incoming particle absorbed in the classical limit 
with probability ^[incoming] = 1. The emission process is described instead by case 1. As 
explained elsewhere [1] this implies that 

lmW(r) = -Im / dr E = (3.40) 



In the end we get, 



am e 



A(r,v)/2 A'(r H {v),v)' 



T = 



4tt 



(3.41) 



confirming the predictions of the Kodama-Hayward theory (|2.4p . We conclude that our 
dynamical black hole is unstable against the emission of spinor particles. It is essential, for 
this result to hold, that the black hole be slowly evolving on the timescale of the wave, since 
otherwise no meaningful notion of a frequency is available. 
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4 MacVittie cosmological black hole 



The preceding calculations showed that the derivatives of the action were strongly mixed by 
the matrix structure of the Dirac equation. Nevertheless things combined so that only the 
radial derivative of the action was really important. Hence we will not study the full Dirac 
equation in the following, but use instead a shortcut. Writing as before 



W = U exp 



n 



(4.1) 



where U is a slowly varying spinor amplitude, from the Dirac equation we get 

$U + fT 1 + m)U = (4.2) 

In the semi-classical limit the second term dominates, so + m) U = 0; thus the matrix 
(i0I + m) must be singular, or equivalently, 

gTd,Jd v I + m 2 = (4.3) 

The action will have a simple pole at the location of the trapping horizon affecting to outgoing 
modes, so for these we also neglect the mass term. Then for an outgoing particle we find 

d r I = -Ffat^dtl , (4.4) 

where 



F(r,t) = y/A s (r){^A s (r) - rH(t)) . (4.5) 

We pick the imaginary part by expanding this function at the horizon along a future null 
direction, using the fact that for two neighbouring events on a null direction in the metric 
(HU, one has t-t = (2H 2 rl)~ 1 (r - r ), where H = H(t Q ). We find the result 

F(r,t)= ^K(ro) ~ r Q H 2 - (r - r„) = n (r r ) (4.6) 

where this time r$ = rjj(to). From this equation we see that d r I has a simple pole at the 
trapping horizon; hence, making use again of Feynman ie-prescription, one finds 

ImJ = -KK{t )- l u{t ), (4.7) 

where u>(t) = dtl is the energy at time t, in complete agreement with the geometric evaluation 
of the previous sections. However, we stress that a full justification of the given shortcut 
really requires the full spinor amplitudes, as we showed above while discussing the BV black 
hole. 



References 

[1] R. Di Criscienzo, M. Nadalini, L. Vanzo, S. Zerbini and G. Zoccatelli, Phys. Lett. B 657 
(2007) 107. 

[2] M. K. Parikh and F. Wilczek, Phys. Rev. Lett. 85, 5042 (2000). 

[3] G. E. Volovik, Pisma Zh. Eksp. Teor. Fiz. 69 (1999) 662 [JETP Lett. 69 (1999) 
705]; A. Calogeracos an d G.E. Volovik, JETP Lett. 69, 281 (1999); G. E. Volovik, 
|arXiv:cond-mat/9902 1 71| 

[4] S. Shankaranarayanan, T. Padmanabhan and K. Srinivasan, Class. Quant. Grav. 19, 
2671 (2002); E. C. Vagenas, Phys. Lett. B 559, 65 (2003); A. J. M. Medved, Phys. Rev. 
D66, 124009 (2002); T. Padmanabhan, Mod. Phys. Lett. A 19, 2637 (2004). 

[5] K. Srinivasan and T. Padmanabhan, Phys. Rev. D 60, 24007 (1999); S. Hemming and 
E. Keski-Vakkuri, Phys. Rev. D 64, 044006 (2001); Q. Q. Jiang, H. L. Li, S. Z. Yang and 
D. Y. Chen, Mod. Phys. Lett. A 22, 891 (2007). 



7 



[6] A. J. M. Medved and E. C. Vagenas, Mod. Phys. Lett. A 20, 2449 (2005); M. Arzano, 
A. J. M. Medved and E. C. Vagenas, JHEP 0509, 037 (2005); R. Banerjee and B. R. Ma- 
jhi, larXiv:0801 .02001 [hep-th] . 

[7] R. Banerjee, B. R. Majhi and S. Samanta. larXiv:0801.3583l [hep-th], 

[8] E. T. Akhmedov, V. Akhmedova and D. Singleton, Phys. Lett. B 642, 124 (2006); 
T. K. Nakamura, arXiv: 0706.2916 [hep-th]; P. Mitra, Phys.Lctt. B648 240 (2007). 

[9] T. Pilling, larXiv:0709.1624 [gr-qc]. 

[10] M. Visser, Int. J. Mod. Phys. D 12 (2003) 649; A. B. Nielsen and M. Visser, Class. 
Quant. Grav. 23 (2006) 4637. 

[11] M. Angheben, M. Nadalini, L. Vanzo and S. Zerbini, JHEP 0505, 014 (2005); M. 
Nadalini, L. Vanzo and S. Zerbini, J. Physics A Math. Gen. 39, 6601 (2006). 

[12] R. Kerner and R. B. Mann, Phys. Rev. D 76, 104010 (2006). 

[13] A. B. Nielsen and J. H. Yoon, larXlv:071 1.14451 [gr-qc], 

[14] S. A. Hayward, Class. Quant. Grav. 15, 3147 (1998). 

[15] H. Kodama, Prog. Theor. Phys. 63, 1217 (1980). 

[16] V. Faraoni, Phys. Rev. D 76 (2007) 104042. 

[17] H. Saida, T. Harada and H. Maeda, Class. Quant. Grav. 24 (2007) 4711. 
[18] R. Kerner and R. B. Mann. larXTv:0710.0612l [hep-th] . 
[19] R. Li and J. R. Ren. larXiv:0802.3954l [gr-qc], 

[20] P. C. Vaidya, Proc. Indian Acad. Sci. A33, 264 (1951); P. C. Vaidya, Nature 171, 260 
(1953); V.V. Narlikar and P.C. Vaidya, Nature 159, 642 (1947). 

[21] J.M. Bardeen, Phys. Rev. Letters 46, 382 (1981). 

[22] C. W. Misner and D. H. Sharp, Phys. Rev. 136, B571 (1964); 

[23] A. B. Nielsen arXiv:0802.3422 [gr-qc]. 

[24] J. R. Ren and R. Li. larXiv:0705.4339l [gr-qc], 

[25] G. C. McVittie, Mon. Not. R. Astronomic Soc. 93, 325 (1933). 

[26] B. C. Nolan, Phys. Rev. D 58, 064006 (1998); Class. Quant. Grav. 16, 1227 (1999); 
Class. Quant. Grav. 16, 3183 (1999). 



8 



